Ab initio four-body calculation of n- 3 He, p- 3 H, and d-d scattering 
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Four-body equations in momentum space are solved for neutron- 3 He, proton- 3 H, and deuteron- 
deuteron scattering; all three reactions are coupled. The Coulomb interaction between the protons 
is included using the screening and renormalization approach as it was recently done for proton- 
deuteron and proton- 3 He scattering. Realistic interactions are used between nucleon pairs. For 
the first time fully converged results for the observables pertaining to the six different elastic and 
transfer reactions are obtained and compared with experimental data. 
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PACS numbers: 21.30.-x, 21.45.+V, 24.70.+S, 25.10.+S 

Modern studies of nuclear structure and reactions re- 
quire ab initio calculations of the quantum many-body 
system using realistic interactions. While at present 
there is a number of successful ab initio structure cal- 
culations using Green's function Monte Carlo (GFMC) 
methods P, [3] and No Core Shell Model ab initio 
studies of nuclear reactions are still limited to few specific 
few-nucleon reactions. The difficulties associated with 
the treatment of the long range Coulomb interaction re- 
stricted, until recently, precise calculations to neutron- 
deuteron (n-d) elastic scattering and breakup Ij], and 
proton-deuteron (p-d) [5j, n- 3 H [fj and p- 3 He [7[ elastic 
scattering at low energy. The situation has now changed 
due to the work in Ref. [8|. Using the ideas of screen- 
ing and renormalization [9( fully converged results for 
observables in p-d elastic scattering and breakup were 
obtained for a wide range of energies and configurations 
using realistic force models. Other developments took 
place recently using the GFMC's approach [lOj to cal- 
culate the n- 4 He phase shifts at energies well below the 
inelastic threshold. 

In the present work we attempt to further extend ab 
initio calculations of the four-nucleon (47V) system by 
solvingthe Alt, Grassberger and Sandhas (AGS) equa- 
tions for the n- 3 He, p- 3 H and d-d reactions. This is 
the most complex scattering calculation that has been 
attempted so far since all three reactions are coupled 
and involve both isospin T = and T = 1 states, to- 
gether with a very small admixture of T = 2 due to the 
charge dependence of the hadronic and electromagnetic 
interactions. This work is the continuation of two pre- 
vious works for n- 3 H [l2| and p- 3 He [l]| where the 4JV 
scattering problem was calculated with the same level of 
accuracy as it already exists for three-nucleon (3iV) sys- 
tem. This means that calculations are carried out with- 
out approximations on the two-nucleon (2N) transition 
matrix (t-matrix) like in Ref. 14| or limitations on the 
choice of basis functions as in Ref. 15]. Therefore, af- 
ter partial wave decomposition, the AGS equations are 
three-variable integral equations that are solved numeri- 
cally without any approximations beyond the usual dis- 
cretization of continuum variables on a finite momen- 



tum mesh. Therefore discrepancies with data may be 
attributed solely to the underlying 2N forces or lack of 
3N forces. 

The equations we solve are based on the symmetrized 
four-body AGS equations of Ref. [13]. In order to in- 
clude the Coulomb interaction we follow the methodology 
of Refs. [H, HH and add to the nuclear pp potential the 
screened Coulomb one wr that, in configuration space, 
is given by 



wr(t) = w(r) e 



(r/Ry 



(1) 



where w(r) — a e /r is the true Coulomb potential, a e ~ 
1/137 is the fine structure constant, and n controls the 
smoothness of the screening; n — 4 to 8 are the op- 
timal values that ensure that Wu(r) approximates well 
w(r) for r < R and simultaneously vanishes rapidly for 
r > R, providing a comparatively fast convergence of the 
partial-wave expansion. The screening radius R must be 
considerably larger than the range of the strong inter- 
action but from the point of view of scattering theory 
wr, is still of short range. Therefore the equations of 
Ref. [13] become R dependent. The transition operators 
where a{(3) — 1 and 2 corresponds to initial/final 
1 + 3 and 2 + 2 two-cluster states, respectively, satisfy the 
symmetrized AGS equations 
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for 2 + 2 as the initial state. In both sets of equations 
Go is the four free particle Green's function and the 
27V t-matrix derived from nuclear potential plus screened 
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Coulomb between pp pairs. The operators U? R ~. obtained 
from 

U" R) = PoGo 1 + P a tW Go Ufa , (4a) 
P\ = P\2 P23 + P13 P23, (4b) 

P-2=Pl 3 P24, (4C) 

are the symmetrized AGS operators for the 1 + (3) and 
(2) + (2) subsystems and Py is the permutation operator 
of particles i and j. Defining the initial/final 1 + (3) and 
(2) + (2) states with relative two-body momentum p 
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the amplitudes for all possible two-cluster transitions are 
obtained as (p/|T ( ^|pi) = S a /3(0a (p/)|W ( ^|<^ fl) (pi)) 

with Sn = 3, S21 = V3, S22 = 2 and S12 = 2\/3. 

In close analogy with p-d elastic scattering, the full 
scattering amplitude, when calculated between initial 
and final p- 3 H or d-d states, may be decomposed as fol- 
lows 
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FIG. 1: (Color online) Differential cross section and neutron 
analyzing power of elastic n- 3 He scattering at 1, 3.5, and 
3.7 MeV neutron lab energy. Results obtained with poten- 
tials CD Bonn (solid curves) , AV18 (dashed curves) , INOY04 
(dashed-dotted curves), and N3LO (dotted curves) are com- 
pared. The cross section data are from Ref. [2l[, A y data are 
from Ref. [13] at 1 MeV and from Ref. [H at 3.7 MeV. 



with the long-range part V^r being the two-body t- 
matrix derived from the screened Coulomb potential of 
the form between the proton and the center of mass 
(cm.) of 3 H for a = 1 and between the cm. of both 
deuterons for a = 2. The remaining is the Coulomb- 
distorted short-range part [T^ — t^R 8 a p] as demon- 
strated in Refs. 0, Applying the renormalization 
procedure, i.e., multiplying both sides of Eq. ^ by the 
renormalization factors [Z R }~^ on the left and 
on the right [fj, in the R — ► oo limit, yields the full 
transition amplitude in the presence of Coulomb 
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where the [Z R ]~ {pf\t^\pi) converges (in general, 
as a distribution) to the exact Coulomb amplitude 
(P/I^qc Ip*) between the proton and the cm. of the 
3 H nucleus (between the cm. of both deuterons) and 
therefore is replaced by it. The renormalization factor 
is employed in the partial-wave dependent form as in 
Refs. 



loved 

I El 



A R 



-2i(<TZ-vZi 



(8) 



with the diverging screened Coulomb p- 3 H (d-d) phase 
shift r]1 R corresponding to standard boundary conditions 
and the proper Coulomb one referring to the logarith- 
mically distorted proper Coulomb boundary conditions. 
Obviously, there is no long-range Coulomb force in the 
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FIG. 2: (Color online) Differential cross section and proton 
analyzing power of elastic p- 3 H scattering at 4.15 MeV proton 
lab energy. Curves as in Fig. [Tj The data are from Ref. [24[. 



3 He states; in that case 
and er£ = rfl R = 0. The second term in Eq. 0, 
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after renormalization by [Z R }~2 [Z R J~z, represents the 
Coulomb-modified nuclear short-range amplitude. It has 
to be calculated numerically, but, due to its short-range 
nature, the R —> oo limit is reached with sufficient ac- 
curacy at finite screening radii R as demonstrated in 
Refs. B El 

for p-d and p- 3 He scattering. As found there, 
one needs larger values of R at lower energies, making the 
convergence of the results more difficult to reach. Never- 
theless, for E p > 1 MeV or E d > 1 MeV the method leads 
to very precise results. Depending on the reaction and 
the energy, we obtain fully converged results for n- 3 He, 
p- 3 H, and d-d observables with R ranging from 10 to 20 
fm. 

The results are also fully converged with respect to the 
partial- wave expansion. The calculations include isospin- 
singlet 2N partial waves with total angular momentum 



3 



1(T 



.o 
E 



C, 

^5 



N3L0 
INOY04 
AV18 
CD Bonn 




40 



50 100 150 
©cm. ( de 9) 



50 100 150 

©cm. (deg) 



FIG. 3: (Color online) Differential cross section and deuteron 
tensor analyzing power T20 of elastic d-d scattering at 3 MeV 
deuteron lab energy. Curves as in Fig. [T] The cross section 
data are from Ref. [2a| and T20 data are from Ref. I26I 1. 



w 80 



§40 



0.4 



0.0 



E p = 2.48 MeV 

N3LO Va- ~ 

INOY04 17/ 

AV18 - -- M? 
; .^CD Bonn Jf' 


E p = 6 MeV j?2 

j • 
/• 
/• 


■ ■ ■ ■ i- 

//V \ 
\ 

f' % \ 
^* \ 
\ 


I ■ i i i I ■ i i i I i i 

• ■ 



50 



100 

. (deg) 



150 



50 



100 

, (deg) 



150 



FIG. 4: (Color online) Differential cross section and proton 
analyzing power of p + 3 H — > n + 3 He reaction at 2.48 and 
6 MeV proton lab energy. Curves as in Fi g. [T] The cross 
section data are from Refs. I27II (circles) and [28j ] (squares) at 
2.48 MeV, and from Ref. [2J at 6 MeV. A y data are from 
Ref. [U at 2.48 MeV and from Ref. [H at 6 MeV. 



I < 4 and isospin-triplet 27V partial waves with orbital 
angular momentum l x < 6, 3iV partial waves with spec- 
tator orbital angular momentum l y < 6 and total angular 
momentum J < 4 s -, 4iV partial waves with 1+3 and 2+2 
orbital angular momentum l z < 6, resulting up to about 
15000 channels for fixed 4iV total angular momentum and 
parity. Initial/final two-cluster states with orbital angu- 
lar momentum L < 5 are included for the calculation 
of observables. For some reactions, e.g., n- 3 He elastic 
scattering, the partial-wave convergence is considerably 
faster, allowing for a reduction in the employed angular 
momentum cutoffs. 

The two-nucleon interactions we use are charge- 
dependent (CD) Bonn 17], AV18 [Hj], inside nonlocal 
outside Yukawa (INOY04) potential by Doleschall [l9j . 
and the one derived from chiral perturbation th eory 
at next-to-next-to-next-to- leading order (N3LO) [20(. 
In Figs. Q] through \5\ we show the results for six dif- 
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FIG. 5: (Color online) Differential cross section and deuteron 
analyzing powers ofd + d^p + 3 H and d + d —* n + 3 He 
reactions at 3 MeV deuteron lab energy. Curves as in Fig. [T] 
The cross section data are from Refs. [32l ] (s quar es) and [33| 
(circles). Analyzing power data are from Ref. [33j for d + d — » 
p + 3 H and from Ref. [|| for d + d -> n + 3 He. 



ferent reactions. The most remarkable findings are: a) 
the excellent agreement with the data for the calculated 
d-d elastic differential cross section; b) the lack of a large 
A y deficiency in n- 3 He a nd o - 3 H unlike what was ob- 
served before in p- 3 He [H[l|; c) the overall description 
of d + d — > n + 3 He and d + d — * p + 3 H data, given its 
complex structure. 

Of all the potentials we use only INOY04 represents 



both n- He and p- H thresholds in the correct position 
(-7.72 MeV and -8.48 MeV, respectively). For this reason 
it is perhaps not surprising to see that this potential gives 
rise to the best description of d + d — > n + 3 He and d + 
d — > p + 3 H data, including the differential cross section. 
Nevertheless, even if the correct position of thresholds is 
an important factor, it is certainly not the only one that 
matters since in n- 3 He elastic scattering and p + 3 H — > 
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Work is in progress to study the effect of three-nucleon 
forces through the inclusion of A degrees of freedom. 
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FIG. 6: (Color online) n- 3 He total cross section as function 
of the neutron lab energy. Curves as in Fig. [1] The data are 
from Ref. H. 



n + 3 He INOY04 gives rise to results that are poorer than 
those obtained with other potentials. 

The largest disagreements with data are observed in a 
very small d-d elastic T20 at 90°, and in the p + 3 H — > 
n+ He proton analyzing power and differential cross sec- 
tion at forward and backward angles at low energies. At 
E p = 2.48 MeV we hnd two sets of data that are incon- 
sistent at backward angles but where this work does not 
help to resolve given the wide range of results obtained 
with the potentials we use. At E p — 6 MeV the results 
are almost independent of the choice of potentials but 
they all miss the data at backward angles. The prob- 
lems associated with n- 3 He and p- 3 H reactions may be 
associated with the + excited state of the alpha particle 
that sits between the two thresholds and whose position 
changes with the potential we choose. 

Finally in Fig. [6] we show the total n- 3 He cross section 
as a function of energy. Again we hnd an agreement with 
data that, depending on the potential, is either similar 
(INOY04) or better (N3LO and AV18) than observed in 



Ti hj. As 



3 H total cross section [12f, the curves 



pertaining to AV18 and N3LO calculations exchange po- 
sition at about the same excitation energy relative to the 
n- 3 H threshold (E ~ 1.3 MeV) but now N3LO leads to 
the highest total cross section at low energy instead of at 
the resonance peak. The reason for this behavior is un- 
clear at this time, but may reside in the relative position 
of the 0~ ( 3 Po) inelastic resonance (3(| that is associated 
with the 0~ T = state at about 0.4 MeV excitation 
above the n- 3 He threshold. 

In conclusion, we show for the first time results of ab 
initio calculations involving four-nucleon reactions initi- 
ated by either n- 3 He, p- 3 H, and d-d. Realistic nuclear po- 
tentials are used between hadron pairs plus the Coulomb 
interaction between protons. Although specific observ- 
ables show discrepancies with data that need further in- 
vestigation, namely the inclusion of three-nucleon forces, 
the results give a reasonable overall description of the 
data for all six reactions that show an improvement rel- 
ative to what is found for T = 1 observables 12, 1 1 311 . 
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